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ABSTRACT 

We use the number density of peaks in the smoothed cosmological density field taken 
from the 2dF Galaxy Redshift Survey to constrain parameters related to the power 
spectrum of mass fluctuations, n (the spectral index), dn/dlnfc (rolling in the spectral 
index), and the neutrino mass, m„. In a companion paper we use N-body simulations 
to study how the peak density responds to changes in the power spectrum, the presence 
of redshift distortions and the relationship between galaxies and dark matter halos. 
In the present paper we make measurements of the peak density from 2dF Galaxy 
Redshift Survey data, for a range of smoothing filter scales from 4—33 h- -1 Mpc. We use 
these measurements to constrain the cosmological parameters, finding n = 1.36+Qg|, 
m v < 1.76eV, dn/dlnk= -0.012±g;^, at the 68 % confidence level, where m„ is the 
total mass of three massive neutrinos. At 95% confidence we find m v < 2.48 eV. These 
measurements represent an alternative way to constrain cosmological parameters to 
the usual direct fits to the galaxy power spectrum, and are expected to be relatively 
insensitive to non-linear clustering evolution and galaxy biasing. 



INTRODUCTION 



The space density of peaks in a cosmological density field 
smoothed with a filter is sensitive to the shape of the lin- 
ear power spectrum of mass fluctuations, even for filter sizes 
where the fluctuations are in the non-linear regime (Croft & 
Gaztahaga 1998). In De & Croft (2007), hereafter Paper I 
we explored the sensitivity of the peak density to parame- 
ters related to the initial power spectrum, as well as redshift 
distortions and variations in the galaxy halo occupation dis- 
tribution (e.g., Berlind & Weinberg 2002). The theory of 
peaks in a Gaussian density field was set out in detail by 
Bardeen et al. (1986, hereafter BBKS). The relationship be- 
tween the peak density and power spectrum in BBKS can 
be used to explore such parameters as the spectral index 
and its dependence on scale (Kosowsky and Turner, 1991) 
along with the neutrino mass. In this paper we use data from 
the 2dF Galaxy Redshift Survey (hereafter 2dFGRS, Colless 
et al. 2001), to constrain cosmological parameters using the 
peak density. 

In a recent analysis of the 2dFGRS final data set Cole 
et al. (2005) employed a direct Fourier method to compute 
the power spectrum. These authors put constraints on sev- 
eral parameters using the directly measured power spectrum 
shape. They assumed a primordial form for P(k) with n = 1 
and h — 0.72 along with a negligible neutrino mass. This 
gave preferred values of Q m h = 0.168 ± 0.016 and a baryon 
fraction of ^ = 0.185 ± 0.046( la errors). This analysis 
therefore implies a significantly lower mass density com- 



pared to the Q m = 0.3 often taken as standard. In com- 
bination with CMB data from WMAP (Spergel et al. 2003) 
Cole et al. find Q m = 0.231 ±0.021 Also on large scales some 
evidence was seen by Cole et al. of the baryon oscillations 
predicted by CDM models. Our present work is complimen- 
tary to this large scale linear theory analysis, with the peak 
density enabling constraints to be placed on the linear power 
spectrum and cosmological parameters from data on smaller 
scales. 

This paper is the second in a series. In paper I, we 
examined the number density of peaks in cosmological sim- 
ulations and compared to results from linear peak theory 
(Bardeen et al. 1986, hereafter BBKS) over a range of den- 
sity field smoothing filter scales. This provided knowledge of 
the length scales for which agreement between theory and 
simulation can be expected. The dark matter simulations 
used to compare to peak theory showed good agreement for 
filter scales between 3-30 h~ 1 Mpc. This was assuming that 
the mean interparticle separation was smaller than the filter 
scale (for agreement at better than the 5% level). In ad- 
dition to simple tests using simulations, we created galaxy 
catalogues using lists of dark matter halos and the Halo Oc- 
cupation distribution formalism (Zheng et al. 2005.) Good 
correspondance was found between the peak density mea- 
sured from the galaxies and dark matter for filter scales 
> 4 /i _1 Mpc. 

The peak density can be used to constrain cosmology 
through its dependence on the power spectrum of mass flue- 
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tuations, P(k). In particular we work with the asymptotic 
number density of maxima, i.e. the number density of peaks 
of all heights and which is found by BBKS to be: 

29 - 6^/6 



(1) 



where _R* is defined to be the following ratio of moments of 
the power spectrum: 



Here 
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The scale dependence of P(k) is probed by smoothing the 
density field with Gaussian filter with comoving radius n, 
i.e. multiplying P(k) by e~ (fcrf)2 . Croft & Gaztanaga (1998) 
found that there is an simple relationship that holds to high 
accuracy between the local slope of the power spectrum n 
at wavenumber k and the filter scale rj ~ OAn/k. The peak 
number density is therefore sensitive to the power spectrum 
slope (but not its amplitude) and the parameters which gov- 
ern it can be tested. 

Paper I focussed on several such aspects of the power 
spectrum such as the neutrino mass, the possible rolling of 
the spectral index n, and redshift distortions. It was seen 
that upon increasing the neutrino mass the peak density 
decreases, as it does when the rolling of the spectral in- 
dex is made more negative. The effect of redshift distortions 
was also quantified, showing that they act to suppress the 
number density of peaks on small scales. This information 
is useful for the present paper which deals with the peak 
density in redshift space measured from an observed galaxy 
catalogue. 

Our approach in this paper is to first test our recov- 
ery of the true peak density using mock catalogues derived 
from simulations that have had the 2dF masks and selection 
function applied to them. We then apply the same estima- 
tion techniques to the 2dF data, using the recovery from 
mock data to estimate the reliability of the method. The 
2dF peak density as a function of filter scale is then used to 
constrain cosmological parameters. 

The plan for the paper is as follows: In Section 2 we 
describe the simulated observations (mock catalogues), in- 
cluding the detailed reasoning behind their use. We give an 
overview of the 2dFGRS and in Appendix A describe how 
the mock catalogues were generated. The latter includes a 
description of the selection function, magnitude limits and 
geometrical limitations of the survey and how they were ap- 
plied to N-body simulations to produce mock catalogues. In 
Section 3 we detail tests of the process of finding peaks in 
the simulations and mock catalogues, including selection of 
the volume limits for a given filter size and how we deal with 
incomplete sky coverage. In Section 4 we describe the cal- 
culation of the observed peak density from the 2dF survey 
data using the completeness functions calibrated from mock 
catalogues. We also describe various sources of uncertainty, 
and how they propagate into our evaluation of cosmological 
parameters. In Section 5 we present the best fit values of 
m v , n and dn/dlnfc including error estimates, using Markov 
Chain Monte Carlo analysis. In Section 5 we conclude and 
discuss possible future work. 



2 SURVEY CATALOGUES 

The main aim of this paper is the extraction of accurate 
information about particular cosmological parameters, the 
neutrino mass, slope of the power spectrum and variation of 
spectral index with scale (rolling.) We can make good esti- 
mates of these quantities if we marginalize over parameters 
whose values are already well known. This involves use of 
present day knowledge regarding certain other parameters 
such as the density of matter, cosmological constant, and 
baryon density, Q m h 2 ,QaIi 2 , Q^h 2 . 

We work with the space density of peaks in the 
smoothed galaxy density field. To recover the actual space 
density of peaks we need to address the issues related to 
the limitations and particular nature of the galaxy survey. 
In the 2dF survey data, limitations arise for several reasons. 
For example certain galaxies may not make it into the survey 
because they are faint, but the magnitude limit varies over 
the survey's angular extent. Very bright objects affect the 
detection of those nearby, and so there are exclusion holes 
in the survey. The survey geometry also has a complicated 
nature, making it different from the uniform contiguous vol- 
ume best suited for finding and quantifying peaks. 

To be able to estimate the space density of peaks in 
the Universe from the observations it is therefore impor- 
tant to be able to compute a completeness function for each 
sub volume of the survey region. We do this by computing 
the ratio of the number of peaks in simulations with full 
sampling and no boundary effects to the peak density esti- 
mated from mock catalogues covering the same simulation 
volume. The mock catalogues are simulated observations 
(simulations combined with all the relevant observational 
constraints) . 

When using the mock catalogues, we make the assump- 
tion that the completeness function which we generate using 
them can be applied to the real Universe. This means that 
we assume that the simulations which are the basis of the 
mock catalogues are a close enough approximation to the 
statistical distribution of matter in the real Universe that 
the completeness is valid. We show in Section 3.1 that in 
practice this is not an important restriction as the complete- 
ness is in fact unity (i.e. no peaks are missed) for most of the 
survey volume we actually use. In Appendix A we describe 
how we generate the mock catalogues. As they are created 
to have the same geometry and limitations as the 2dFGRS 
we first describe the observational dataset. 



2.1 Description of the survey 

The 2dFGRS (Colless et al. 2001) resulted in one of the 
largest catalogues of galaxy redshifts made so far, follow- 
ing a major spectroscopic survey which took advantage of 
the unique capabilities of the 2dF facility at the Anglo- 
Australian Observatory. The 2dFGRS obtained spectra for 
245591 objects, mainly galaxies, brighter than a nominal 
extinction-corrected magnitude limit of bj = 19.45. Re- 
liable (2dF quality flag > 3) redshifts were obtained for 
221414 galaxies. The galaxies cover an area of approximately 
1500 square degrees selected from an extended version of the 
APM Galaxy Survey (Maddox et al. 1996) in three regions: 
an NGP strip, an SGP strip and random fields scattered 
around the SGP strip. In figure[T]we show a projection of the 
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Figure 1. 2D projection onto the Cartesian x — y equatorial plane of the 2dF Galaxy Redshift Survey data. The points indicate 2dF 
galaxy positions, and the circles represent peaks in the galaxy density field when smoothed with a 20 h _1 Mpc radius filter. 



galaxy distribution in the 2dFGRS in the equatorial plane 
as well as the positions of peaks found after smoothing the 
density field with a Gaussian filter of radius 20 ft -1 Mpc. 
(the peak finding will be described later). 

Different masks are used when interpreting the 2dF 
data to characterize the completeness of the survey in var- 
ious ways as a function of position on sky. The magnitude 
limit mask gives the extinction-corrected magnitude limit 
of the survey at each position on sky. The redshift com- 
pleteness mask quantifies the fraction of galaxies above the 
magnitude limit with measured redshifts. We use the rec- 
ommended method (Colless et al. 2001) to define this lat- 
ter mask, it being specified by the complete set of 2 degree 
fields that were used to tile the survey region for spectro- 
scopic observations. Each point on the sky inside the survey 
boundary is covered by at least one 2 degree field, but more 



often by several overlapping fields. A sector is defined as 
the region delimited by a unique set of overlapping 2 degree 
fields. This is the most natural way of partitioning the sky, 
as it takes into account the geometry imposed by the pat- 
tern of 2 degree fields and the way in which the galaxies were 
targeted for spectroscopic observations. Within each sector 
with position angle 9,5, the redshift completeness R(8,5), 
is the ratio of the number of galaxies for which redshifts 
have been obtained, N z (6,8), to the total number of ob- 
jects contained in the parent catalogues, N P (8,S) so that 
R(9,S) = N Z {9,5)/N P (9,5). 

We note that when analyzing the data, the redshift com- 
pleteness of a given sector, R(8,S), should be clearly dis- 
tinguished from the redshift completeness of a given field 
denoted, by Cf, since multiple overlapping fields can con- 
tribute to a single sector. The ^t-mask gives the dependence 
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of the redshift completeness on apparent magnitude. We de- 
scribe the /i-mask in more detail in the Appendix A. We use 
these masks (i.e. values of completeness for different patches 
of sky) along with the incorporation of redshift limits, mag- 
nitude limits and a completeness cut-off to make our mock 
catalogues from simulations (see Appendix A). 



3 PEAKS IN MOCK CATALOGUES AND 
SIMULATIONS 

In order to get a feel for how well peaks that are seen in 
the fully sampled dark matter density field can be recov- 
ered in the mock galaxy catalogues, we carry out a one to 
one comparison. To find peaks operationally, we first assign 
the galaxy (or particle) positions to a grid, which is then 
smoothed with a Gaussian filter. Peaks are then associated 
with local maxima in this smoothed density field (grid cells 
which have a higher density than the surrounding 26 cells). 

The first step towards comparing peak locations in 
mocks and the fully sampled dark matter density field is 
to ensure that the same piece of the simulation is used in 
both the mock catalogues and the fully sampled case. As 
mentioned before, our simulations have a box size of 300/i" 1 
Mpc and we replicate this box, making use of the periodic 
boundaries to get the desired survey volume for the mock 
catalogues. We apply the 2dF mask and selection function, 
resulting in a volume that matches the shape of the 2df red- 
shift survey. 

When finding peaks from the mock catalogues it is im- 
portant to assign weights to the galaxies, the appropriate 
ones being related to the inverse of the selection function. 
In order to find the weights, we generate random catalogues, 
in which points are initially randomly distributed. We apply 
the 2dF survey mask and selection function to these points, 
in the same way as with the mock galaxy catalogues. We 
then place a grid on top of this distribution of points and 
count the number of points in each grid cell. The ratio of the 
number of points in the grid cell to the expected number for 
a uniform distribution gives the inverse of the appropriate 
weight to use for that cell. In this way, we compute a weight 
for each cell on the grid. In order to make sure the weights 
have an acceptably low contribution from Poisson noise we 
make the number of points in the random catalogues twenty 
times that in the mock galaxy catalogues. 

We place the same three-dimensional grid (of size 256 3 
cells) onto the distribution of galaxies in the mock cata- 
logues. The box size was set equal to the distance R ma x as 
described earlier. We assign the weighted galaxy density to 
the grid and smooth this field in Fourier space with a Gaus- 
sian filter. We carry this out for several different smoothing 
filter radii, r/, each time locating the local maxima in the 
field. 

For our one-to-one comparison, we compare directly the 
positions of peaks found from the mock catalogues and from 
the fully sampled simulations. Because of shot noise, we 
don't expect peaks in the mocks to be found in exactly the 
same places. In Figure [2] we plot projections of the NGP 
and SGP regions for one of the mock catalogues, with sym- 
bols showing where the peaks were found and other sym- 
bols showing where the true peaks lie (measured from the 
full simulation.) We can see there that the distrbutions of 
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Figure 2. A test of peak finding in simulated data. We show a 2D 
projection of the positions of peaks found from a fully sampled 
simulation (circles) and a mock catalogue (crosses) for a filter 
scale Tf = 10.2/i — 1 Mpc. The top panel shows the SGP and the 
bottom panel the NGP. 
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Figure 3. Percentage of peaks found from mock catalogues (of 
the NGP region) within 1 smoothing filter radius of peaks mea- 
sured from the fully sampled simulation. We also show the re- 
sults for peaks measured from random catalogues. Both results 
are shown as a function of filter radius r * . 
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points are similar in general, but that there are differences, 
particularity at the far edges of the survey. 

We quantify this in Figure [3] where we show the frac- 
tion of peaks in the mock catalogues which have a true peak 
within one filter radius. We can see that good peak detec- 
tions, defined in this way occur approximately 50% of the 
time, roughly independent of the filter size. In the same plot 
we have carried out this comparison with peaks found from 
the random catalogues, to show how likely it is that false 
peaks arise from Poisson fluctuations. We can see that this 
is ~ 10 — 20 times less likely, showing that we are truly 
extracting peak information from the mock catalogues. 

The one-to-one test carried out in this section is meant 
to be illustrative of how well peaks can be identified in 
galaxy catalogues. The information presented in Figures [5] 
and [3] is not used in our measurement of the peak density 
from observations. The completeness corrections used for 
this purpose are described in the next section. 



3.1 Determination of an effective peak 
completeness 

In order to compare an observed peak density with theoret- 
ical models, we need to statistically correct for the effects of 
shot noise and of survey boundaries. Our one-to-one com- 
parison above does not show the aggregate effect on the peak 
density of these effects, but only that approximately 50% of 
peaks are found within 1 filter radius of their true locations. 
If we make the assumption that the 2dF observations are af- 
fected in the same way by boundaries and shot noise as our 
mock catalogues, we can use our mock catalogues to com- 
pute an effective completeness (which varies spatially) . This 
will enable us to translate an observed peak density into the 
density we would have measured with fully information and 
denser sampling. 

Our calculations of the effective completeness involve 
first locating the local maxima in the simulation with full 
sampling and in the mock catalogues, as described above. By 
examining the peak density as a function of distance from 
the edge to the survey we have found that boundaries do 
have a substantial effect (for example, in Figure [T] it can be 
seen that peaks occur preferentially at the edges) . Rather 
than applying a large correction for this effect, we decide to 
not consider the peaks in the volume near the boundaries. 
We choose to remove from our number density computation 
all peaks and volume closer to the side edges of each bin 
within 1.5rf/R ma x in azimuthal distance and rf/R max in 
declination distance on each side. Our approach is to only 
include regions of the simulation for which the completeness 
correction is of order unity (see below.) 

We compare the number of peaks from fully sampled 
simulations and mock catalogues to yield the completeness. 
To do this we divide space into radial bins (distance mea- 
sured from the origin) of width 0.5r/ for filter radii less than 
15/i~ Mpc and of width of O.lr/ for larger filter radii. For 
each radial bin i, a completeness d given by 
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Figure 4. The peak density measured from the NGP (bottom 
panel) and SGP (top panel) mock catalogues. The long-dashed 
lines show the average for the mocks (we plot the corrected peak 
density when showing the mock results.) The points are results for 
each of the 20 mock catalogues. The short-dashed lines indicate 
the peak densities measured from the fully sampled simulations 



(nstm)i and (n moc k)i are respectively the number of peaks in 
the i th bin measured from the simulation and the mock cat- 
alogues. The completeness correction for a particular filter 
size, Tj (denoted by Cj(r/)) is found by averaging over com- 
pleteness values obtained from 20 mock catalogues. From an 
observed number of peaks, the corrected number of peaks in 
the i th bin, rij can therefore be computed from 



H = Ct{rf) x n obs 



(4) 



where n b s ,i is observed number of peaks in the i* radial 
bin for the same filter radius. 

For most of the mock survey volume, the completeness 
we which we compute in this way is Cj = 1.0, indicating 
that we are statistically finding all peaks, with no need for 
any correction. The completeness correction therefore has 
a small effect, and there is no significant difference in our 
results if we do not include it. The completeness correction 
does however enable us to include a larger volume of the sur- 
vey than without it, resulting in smaller error bars. The com- 
pleteness correction d{rj) rises above 1.0 for large distances 
from the origin, as the number density of galaxies goes down. 
For example, with a filter scale 77 = 8 ft -1 Mpc between a 
radial distance of 300 h~ Mpc (where C\(rf) ~ 1.0) and 
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Figure 5. A comparison of the redshift distribution of galaxies in 
one of the mock catalogues and the 2dFGRS. We show results for 
the NGP and SGP regions separately. The version of the mocks 
and 2dFGRS used here was that appropriate for a smoothing 
filter scale rf = 9.5 h~ Mpc (i.e. the RA and dec cuts were those 
appropriate for that rf, as descibed in Section 3.3 of the text.) 



the edge of the survey volume used (at 360 h~ Mpc) Ci(r/) 
averages approximately 1.7. 

In Figure [3] we have plotted peak densities for differ- 
ent 77 values for both NGP and SGP mock catalogues. We 
have also shown with lines the average for mock catalogues. 
The averaged corrected peak densities (multiplied by r^) are 
shown for both the NGP and SGP mocks. We find that for 
all r-f values they agree very well with the peak densities 
from the fully sampled simulations. 

In Figure [5] we have shown the number of galaxies in 30 
radial bins for r/ = 9.5/i _1 Mpc. We have already seen that 
in the very nearby universe it was necessary to superimpose 
two mock catalogues in order to achieve the required number 
density. 
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Figure 6. The peak number density measured from 2dF obser- 
vations as a function of filter scale rf. We show results for the 
NGP and SGP regions separately as symbols. We also show a 
line correponding to the best peak theory fit (see Section 5 for 
details.) 
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Figure 7. The observed peak density measured from the 2dFGRS 
(points) alongside several different peak theory model curves. The 
2dFGRS results are the weighted average of those for the NGP 
and SGP regions show in Figure [6] The models are examples of 
number densities computed from peak theory for a few different 
values of the spectral index, n, the running of the index, d and 
the total mass of neutrinos, m v . 



4 PEAK DENSITY FROM THE 2DFGRS 

We use Equation [4] to calculate the peak density in the ob- 
served universe from the 2dFGRS data. To compute the 
correction factors, we use 20 mock catalogues which have 
a comparable number density of objects to that in the sur- 
vey. These mocks are used to calculate the error bar on each 



data point corresponding to each rf. In Figure|6]we show the 
peak densities for the SGP and NGP. We also indicate the 
best-fit theoretical curve which corresponds to a model with 
n = 0.94, d = -dn/dlnk = -0.009, m„ = 0.70 eV (total 
mass of three massive neutrino species in eV) . We descibe 
how the model fitting was done in Section 5 below. 

In Figure [7] we show the weighted average of the peak 
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Figure 8. The covariance matrix of the peak density com- 
puted using mock catalogs and Gaussian simulations (see §5.2)). 
The symbol area is proportional to dj /(CiiCjj) 1 / 2 , with nega- 
tive elements shown as open symbols. The 30 elements span values 
of r-f = AA7 ?i~ 1 Mpc to r f = 33.1 h _1 Mpc, as given in Tablc[l] 



density measured from the SGP and NGP. For illustrative 
purposes, we also show the peak densities predicted by vari- 
ous models. The observational weighted average of NGP and 
SGP data was calculated using: 



^pk.obs 



(/sGp)^pfc,NGP + (/NGp)^pk,SGP 



(5) 



/NGP + /SGP 

where n p k,sGP and /sgp are respectively the number density 
of peaks in the SGP region of the survey and the fractional 
error on that value. 

In order to carry out the maximum likelihood fitting de- 
scribed in the next section we compute the error in our mea- 
surement of the peak density from observations. We make 
use of our mock catalogs to quantify the different compo- 
nents. 

We include the contribution caused by cosmic variance. 
This we estimate from the standard deviation about the 
mean in the number of peaks in the different mock cat- 
alogues. We compute this for the 25 different filter radii 
ranging from 4.5 /i _1 Mpc to 33 /i _1 Mpc, using the 20 dif- 
ferent mock catalogues. We also use the scatter between the 
mock catalogues to quantify the effects of variations in the 
completeness corrections on the final result. For the error 
associated with the weighted average of the NGP and SGP 
regions, we use error propagation, finding the result to be 
very close to |(cr^ 9p + a 2 gp )2 . 



5 CONSTRAINTS ON COSMOLOGICAL 
PARAMETERS 

5.1 CDM parameters 

In order to see which values of cosmological parameters are 
allowed by the data, we fit the peak density we measure 
from the 2dFGRS observations with that expected in dif- 
ferent versions of the CDM cosmology. We concentrate on 
three parameters, the spectral index n, the running of n, 
dn/dlnfc, and m v (the total mass of three massive neutrinos 
species in eV). For the other parameters which affect the 



Table 1. The space density of peaks measured from the 2dF 
redshift survey data, as a function of filter scale rf 



T'pkT f 



(ft _1 Mpc) xlOOO 



o"(n pfc r3) 
xlOOO 



4.47 


4.97 


0.11 


4.79 


5.21 


0.10 


5.13 


5.45 


0.11 


5.50 


5.47 


0.11 


5.89 


5.56 


0.11 


6.31 


5.81 


0.11 


6.76 


5.98 


0.11 


7.24 


6.00 


0.11 


7.76 


6.43 


0.12 


8.32 


6.28 


0.12 


8.91 


6.33 


0.12 


9.55 


6.80 


0.13 


10.2 


6.97 


0.13 


11.0 


7.39 


0.14 


11.7 


7.33 


0.14 


12.6 


7.61 


0.15 


13.5 


7.85 


0.15 


14.5 


7.90 


0.16 


15.5 


7.81 


0.16 


16.6 


8.16 


0.17 


17.8 


8.34 


0.17 


19.1 


8.28 


0.18 


20.4 


8.78 


0.20 


21.9 


9.47 


0.21 


23.4 


8.90 


0.22 


25.1 


8.31 


0.19 


26.9 


9.24 


0.22 


28.8 


9.30 


0.24 


30.9 


8.90 


0.24 


33.1 


9.34 


0.27 



shape of the power spectrum, we use values from the WMAP 
first year data release (Spergel et al. 2003). These are 
n m h 2 = 0.135±g;gg|, h = 0.71±£:o3, ^6 = 0.0224 ± 0.0009. 
Our results are insensitive to adoption of the WMAP5 pa- 
rameters instead (Dunkley et al. 2009) . We have seen in Pa- 
per I that the effect of small scale redshift distortions can be 
parametrized using the one dimensional velocity dispersion 
of galaxies, aj . We use a prior on its value in our analysis: 
07 = 3.4 ± 1.0 (in units of 100 km s _1 ), based on the value 
given Matsubara et al. (2004). The Gaussian error on the 
value we use is large enough to incorporates other recent 
measurements at the la level such as that of Peacock et al. 
(2001) (400 km s _1 ). We also use a prior on the running of 
the spectral index, = ± 0.3. 



5.2 Construction of the Covariance Matrix 

In our analysis we compute the likelihood as oc e~ x ^ 2 , where 
the x 2 values use the full covariance matrix to include the 
correlations between bins at different scales. We note that 
we have found that including only the diagonal elements 
would give erroneously tighter constraints (at the factor of 
~ 2 level) on cosmological parameters. We would like to be 
able to calculate the covariance matrix using our 20 mock 
catalogues, but with only 20 of them it is too noisy to invert. 
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We therefore use the mock catalogs to compute the diagonal 
elements only. We then use 1000 Gaussian realizations of the 
density field to compute an alternative, less noisy version of 
the covariance matrix. We scale the off diagonal elements of 
this matrix so that they are consistent with those from the 
mocks (see below). 

Operationally, each element of the matrix is defined 

from 



' i 



E 



(6) 



where Xi is the value of n p k for 77 bin i and the mean Xi 
is computed from the 20 mock datasets or 1000 Gaussian 
realizations. 

We compute the x 2 of the difference between observa- 
tional measurements and a theory curve using 

X ^ ^ (^i,obs theory ) (C )ij (^j.obs theory) (7) 

In equation [7] 2i, i, s is r-f bin i from the observations and 
similarly Zi,theory indicates the i-th theory bin. The sum is 
over all the r/ bins in the dataset (in our case there are 
30 bins). In using the covariance matrix the main challenge 
is the inversion of dj to calculate % 2 ■ We have chosen the 
singular value decomposition technique to do this. We have 
plotted the elements of dj in graphical form in Figure[8] We 
can see that the non-zero elements are concentrated close to 
the diagonal, as expected, but that there is still significant 
covariance between bins which must be accounted for, so 



mock 
ij 



that we need to invert the matrix and use equation [7] Even 
with 20 mock catalogues, dj is quite noisy, making inversion 
difficult. 

To deal with this, as stated above, we make Gaus- 
sian linear theory simulations (with the same linear the- 
ory power spectrum as our N-body simulations, and the 
same box size). We compute the peak number density 
in these simulations and use them to compute the off- 
diagonal elements in the covariance matrix, using dj 

//nmock /imock \ 1/2 ^Gauss. / /^Gauss. pGauss. \ 1/2 tt ov „ f~t, 

\yu ^jj ) ^ij I v^u ) ■ riere o. 

are elements computed from the mocks and Cjj auss ' from 
the Gaussian simulations. We do this for each half of the 
survey. The resulting matrix elements for the NGP are then 
summed with those for the SGP calculated in a similar fash- 
ion and this matrix is the one used to compute equation [Jj 
in our analysis . 



5.3 Markov Chain Monte Carlo 

In order to compute constraints on the multi-dimensional 
space of parameters efficiently, we use use a Monte Carlo 
Markov Chain approach. The procedure that we follow is 
that of Verde et al. (2003). As in that paper we generate 
Markov chains which have a large number of values (30000) 
to ensure a satisfactory level of convergence. 

To generate the Markov chain we compute theo- 
retical peak densities for an initial parameter set of 



Peaks in the cosmological density 9 




2000 

1800 
1600 
1400 
1200 
1000 
800 
600 
400 
200 




2000 
800 
600 
400 
200 
000 
800 
600 
400 
200 




-0.5 -0.4 



-0.2 -0.1 
dn/dlnk 



0.1 0.2 0.3 




Figure 10. One dimensional probability distribution of cosmological parameters measured from the 2dFGRS data. We show results for 
spectral index n, the running of the index dn/dlnfc, and the neutrino mass m v (in eV). The y-axis on each panel gives the number of 
Monte Carlo Markov Chain points in each bin. 



n, dn/dlnfc, 07, f2„ values chosen randomly (our results are 
not sensitive to the exact range they are drawn from.) We 
compute the x 2 of the fit to the peak density data from 
the 2dFGRS and then the likelihood, £ = e" x2/2 . We then 
randomly pick another parameter set using a Gaussian dis- 
tribution centered on the previous set (i.e. mean and vari- 
ance equivalent to the a 2 associated with each parameter.) 
If the likelihood found from the new set of parameters is 
higher than before, we include this new parameter set in the 
Markov chain and move on, otherwise we compare a ran- 
domly generated number between and 1 to the ratio of 
the new and old likelihoods. If the number is lower than the 
ratio we decide to include the new set of parameters in the 
chain and if it is not we return to the old set of parameters 
instead. As in Verde et al. we start the Markov Chain from 4 
different random locations in parameter space and combine 
the results. 

In figure [§] we show the two dimensional probability 
density contours for the 1,2,3- a (68.3%, 95.4% and 99.7% 
levels.) We can see that there are some strong degeneracies 
between the parameters, with for example the effect of a 
large n being balanced by a larger m„ . Promisingly, the value 
of 07, parametrizing the suppression of peak density due 
redshift distortions does not have a strong effect on the other 
parameters. Overall, we can see that the points with n = 1, 
dn/dlnk — and m„ = lie inside the 1-2 a contours, 
indicating that standard assumptions for these parameters 
are reasonably consistent with the data. 

In Figure [10] we show the 1 dimensional marginalized 



Table 2. The best fit values for three cosmological parameters 
measured from the peak density in the 2DFGRS data (see §5.2) 
for details. 



parameter best 1 a 2 a 

fit range range 



n 1.36 0.72^2.12 0.14 -* 2.61 

dn/dlnk -0.01 -0.22 -> 0.18 -0.33 -» 0.20 
m„ (eV) 0.27 0.0 -» 1.76 0.0 -> 2.48 



probability distributions for the 3 parameters n, dn/dlnfc, 
m„. The width of the likelihood distributions for n, dn/dlnk 
are large, indicating that the peak density constrainty will 
not be competitive with other measures of these parame- 
ters. The m„ probability distributions are relatively narrow, 
however. 

The best fit values of the parameters in the marginal- 
ized distributions and their range of uncertainty are give in 
Table [5] We can see that n is consistent within the large 1 a 
error bars with the WMAP5 results (Dunkley et al. 2009), 
n — 0.963^qqJ^. The 2 a upper limit on the total mass of 3 
neutrino species, m„ is 2.48 eV. This is in the same range 
as upper limits from a variety of other cosmological probes 
(see Elgaroy 2007 for a recent review.) 



10 S. De and R.A.C. Croft 



6 CONCLUSIONS 

By measuring the space density of peaks in the smoothed 
2dF galaxy survey density field as a function of filter scale, 
we have constrained the values of the cosmological power 
spectrum parameters dn/dlnfc, n and Q. v . We find values 
that are consistent within 1 a of those measured from other 
cosmological observables (e.g., Dunkley et al. 2009, Seljak 
et al. 2005), and are in support of the standard ACDM 
scenario. In linear theory, the peak density is not affected 
by the amplitude of fluctuations. As simulation tests have 
shown that this is also true in the non-linear regime, for 
filter scales as small as 3 /i _1 Mpc, this means that our con- 
straints on cosmological parameters come entirely from the 
shape of the power spectrum. Because of degeneracies in the 
shape of the power specrum with cosmological parameters 
it was necessary to assume prior values for some parame- 
ters (fi m , tt b and H ) from the WMAP (Spergel et al. 2003) 
results in order to derive our own limits. 

Various sources of systematic errors are possible, when 
measuring peak statisics, for example those related to the 
fact that the peak density in a smoothed field is sensitive to 
boundary effects, which are difficult to correct for (unlike the 
case of the correlation function, for example.) In the present 
work, we have made conservative cuts, ignoring space close 
to the survey boundaries and validated these using mock 
catalogs to test both the overall recovery of the peak density 
and the recovery of individual peaks. An additional source 
of error comes from the fact that redshift distortions tend to 
merge peaks together which are separate in real space, as was 
seen in paper I. We have dealt with this effect, by including 
a power spectrum suppression term due to the small scale 
random velocity dispersion (see paper I for details and tests) . 
This is the one part of the analysis for which peak finding 
may be most sensitive to the relationship between mass and 
light. Our tests using the halo occupation distribution in 
Paper I have show that this is unlikely to be a problem at 
the level of current observational uncertainties. 

Our constraints have come from an analysis of filter 
lengthscales between 4 — 30 ft _1 Mpc. This approach is 
therefore complementary to much work in large scale struc- 
ture which has been done by analyzing the galaxy power 
spectrum on somewhat larger length scales (wavenumber 
fc < 0.1 ft Mpc -1 , corresponding to wavelengths ~ ir/k > 
30 /i _1 Mpc). On these larger scales, galaxy fluctuations are 
assumed to be linearly related to mass, and direct measure- 
ments of the galaxy power spectrum are used to constrain 
the mass power spectrum and hence cosmological parame- 
ters (e.g., Tegmark et al. 2004). 

Recently, it was shown by Sanchez & Cole (2007) that 
shapes of the power spectra measured from the 2dF and 
SDSS galaxy surveys (e.g., Cole et al. 2005, Tegmark et al. 
2004) are not in agreement, due to the r-band selected SDSS 
galaxies having a stronger scale-dependent bias. This results 
in differences in the cosmological parameters inferred from 
the two surveys which are larger than the quoted measure- 
ment uncertainties. From the tests in CG97 and Paper I, we 
expect that our peak based measurement to be more robust 
both to differences in galaxy bias and non-linear evolution. 
In the future it will be useful to compute the peak number 
density from the SDSS survey data, in order to check consis- 
tency. The larger size of the SDSS final dataset will enable 



the use of larger filter scales and reduce the error bars on 
parameters. 
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APPENDIX A: MOCK CATALOGS 

Al Construction of Mock catalogues 

In order to incorporate observational constraints into a sim- 
ulation we need to first compute the the selection function. 
We choose our luminosity function to have the form of a 
Schechter (1976) function (following Erdogdu et al. 2004) 
and define 



i/>(r) = 



$(L)dL 



$(L)dL 



(8) 



where ip{r) is the radial selection function at a distance 
r which is obtained by integrating over the luminosity func- 
tion $(-£/)• L m i u (r) is the minimum luminosity which could 
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be observed at a distance r. Lo is the lowest measured lumi- 
nosity in the survey, which is taken as the luminosity L m i n (r) 
at the comoving distance r of 5/i _1 Mpc, equivalent to a 
lower redshift cutoff for the data sample we use. 

To choose magnitudes for the galaxies in our mock cat- 
alogues, we numerically solve equation (1) for L min (r) using 
a random value of the selection function between to 1. Us- 
ing this value of the luminosity we estimate the apparent 
magnitude at point (r,9,S), 9 being the azimuthal angle, 8 
being the declination and r being the comoving radial dis- 
tance. There is a dependence on angular coordinates that 
comes from the fact that Lo is estimated as the minimum 
luminosity at a fixed radial distance but its value varies with 
9 and 5. Once we have found the minimum luminosity at r 
from the last integral equation, we use following equation to 
compute the apparent magnitude (m): 

Lmi„(r) = [d L (l + z)] 2 (l + z)W W.O-OA (m -M,) (g) 

where <1l is the luminosity distance, z is the redshift, 
and M* = —19.7, equivalent to the absolute bj limit. The 
extra factor of (1 + z) is due to the inclusion of the k cor- 
rection. 

The mock catalogues are made by applying the mask 
and luminosity function mentioned above to the outputs of 
N-body simulations. To generate the simulation we use the 
same cosmology as described in paper I, an LCDM run- 
ning spectral model with spectral index n = 1, dn/dlnfc = 
—0.020, Q m = 0.3, S7a = 0.7, and without massive neutri- 
nos. We use the transfer function described in Eisenstein and 
Hu (1999). These simulations are identical to those used in 
paper I with respect to box size (300fa -1 Mpc), mass fluc- 
tuation amplitude, and other relevant parameters (particle 
number = 3.375 x 10 6 ). 

In paper I, we tested the effect on the space density 
of peaks of using a galaxy Halo Occupation Distribution 
(HOD, e.g., Berlind & Weinberg 2002) to create the galaxy 
density field rather than particles. We found no measurable 
difference for filter scales 77 4 /i _1 Mpc and therefore in 
the present paper, in the interests of having the maximal 
space density of objects, we use particles as proxies for galax- 
ies in our mock catalogues. 

In order to produce mock catalogues, we first apply the 
geometrical constraints (angular coverage) of the survey to 
our simulations. There are galaxies in the 2dF survey at red- 
shifts as deep as z — 0.3. To cover the whole of this observed 
volume it is necessary to replicate the particle distribution of 
the simulation (box size of 300/i 1_ Mpc) so that it can cover 
a sphere of radius ~ 8OO/1 1 Mpc. To do this, we translate 
all particles in the same direction with a comoving distance 
equal to the box size (the simulations have periodic bound- 
ary conditions), repeating this in three mutually orthogonal 
directions to yield enough volume to inscribe the sphere in 
it. We note that although the box has been replicated several 
times, the actual volume of a single box is similar to that of 
the survey, because of the wedge like geometry of the survey. 
Having carried out this replication, we compute the comov- 
ing distance from the observer for our LCDM cosmology, 
and compute the redshifts of all the particles, including the 
effect of peculiar velocities. 

After we select the simulation particles which have the 
same angular coverage as in the real survey, we reject those 
which have a sector completeness, Cf less than 0.2. We as- 



sign random values to the selection function to calculate 
L m in(r) as described above. This equation is then solved 
to obtain apparent magnitudes for each galaxy. We apply 
magnitude limits (described below) to exclude faint objects 
from the mock catalogues. One important point to note is 
that in the survey there is also a bright magnitude limit. It is 
necessary to ignore very bright galaxies as they can affect ob- 
servations of neighboring galaxies. The presence of holes in 
the angular mask to exclude them also results in holes in our 
mock catalogues. Our chosen values for the bright and faint 
magnitude limits are bj = —15.0 and bj — —19.4 respec- 
tively. Galaxies falling outside this range are ignored. Using 
the computed apparent magnitude we calculate the magni- 
tude dependent incompleteness c z (m,/i), which is given by 

c(m,Ai) = 7 [l-e (m -^] (10) 

where 7 = 0.99 and fi (different for each mask cell) are 
parameters which were set by Colless (2001) by comparing 
to a simple power law model for galaxy number counts (see 
equations 5 — 8 of Colless 2001.) These values of fi were ob- 
tained from masks made publically available by the 2dFGRS 
collaboration For each mock galaxy we compare this com- 
pleteness to a randomly generated number between and 1, 
rejecting galaxies when the number is greater than the com- 
pleteness. After this, the final mock catalogue is output. To 
be able to average over random fluctuations we make many 
catalogues from 5 simulations generated with different ran- 
dom seeds as well as by randomizing the observer positions 
in the simulations. We use 4 different observer positions per 
simulation, making 20 different mock catalogues. 

One problem we encounter when making the mocks is 
that for small distances from the observer the number den- 
sity of particles in the simulations is less than the number 
density of galaxies in the survey. For this region, we clone 
particles to make up the difference. Our results are insensi- 
tive to whether this is done, as we restrict our analysis to 
smoothing filter radii much larger than the particle mean 
separation at this distance. We describe this procedure in 
more detail below. 



A2 Determination of the radial boundary 

Before determining the peak density from the catalogues 
it is necessary to choose the region which can be analysed 
for each given filter size. This is because the space density 
of galaxies decreases with distance from the origin. If the 
space density of galaxies is too low for a given filter size, 
the smoothed density field from the catalogue will not give 
a correct representation of the smoothed underlying density 
field. In paper I, it was shown that the smoothing filter scale 
should be greater than the mean separation between galaxies 
to avoid serious underestimation of the peak density. 

Additionally, we have a finite sized (256 3 cells) grid 
which we use to locate peaks (local maxima on the grid). 
We make sure at all times that the grid cell size is also 
less than the mean separation between galaxies, in order to 
avoid missing peaks. If w is the cell width and r g is the mean 
intergalaxy distance, then we ensure that r fl ^ w. We have 
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Figure 11. Variation of the mean intergalaxy separation r p in 
radial bins centered on distance Rmax from the origin in mock 
catalogues and 2dFGRS observational data (for the NGP half of 
the survey) 

checked to see that this criterion is sufficient, and our results 
are not sensitive to this exact choice. 

For a given value of -Rmax, we construct a volume lim- 
ited catalogue, and compute the mean separation between 
galaxies. In Figure [11] we show the dependence of this aver- 
age inter-galaxy separation r g (~ [V/Ng]?, where V is the 
volume in a radial bin centered on Rmax and N g is the num- 
ber of galaxies) with distance from the origin Rmax in both 
mock catalogues and the actual 2dF survey. 

As mentioned above, for each filter radius r/ we can 
safely increase Rmax to the maximum value where r/ r g . 
In order to determine this value, we averaged over the mock 
catalogues and fitted a fourth order polynomial in log(r/) 
to yield a fitting function for the radial distance, J? ma:E (77). 





o 

CO 



• • #0# • • • • 

•••••• #Q# • • • • 

• • • • • • • • • 

• • • • • 



. • • • 
• • • • • • • 

• • • • • • • • 

• • • • • • • • 

• • • • • • • • • 



o 

C\2 



• • • • • 

• • • • • • 

•• • • • • ' 

•••• #£# • • • • * • 

••• #0# ••••••• 

• • • #0# • • 

• • • • ### • • • 

• • • • • • 

- • • • • • 

•••••• • • 



• • • • • 

• • • • 

• • 

• • 





o 







01 



02 



OG 



